January 4, 2012 2:43 WSPC/INSTRUCTION FILE sulaiman 



International Journal of Modern Physics: Conference Series 
© World Scientific Publishing Company 



MODELING THE INTERACTIONS OF BIOMATTER AND 

BIOFLUID 



A. SULAIMAN 

Badan Pengkajian dan Penerapan Teknologi, BPPT Bid. II (19 th floor), 31. M.H. Thamrin 8, 

Jakarta 10340, Indonesia 
asulaiman@webmail. bppt.go.id, sulaiman@teori.fisika. lipi.go.id 

L.T. HANDOKO 

Group for Theoretical and Computational Physics, Research Center for Physics, Indonesian 
Institute of Sciences, Kompleks Puspiptek Serpong, Tangerang, Indonesia 
handoko@teori.fisika.lipi. go. id, handoko@fisika.ui. acid, laksana.tri.handoko@lipi.go.id 

Department of Physics, University of Indonesia, Kampus UI Depok, Depok 164-24, Indonesia 

Received Day Month Year 
Revised Day Month Year 

The internal motions of biomatter immersed in biofluid are investigated. The interactions 
between the fragments of biomatter and its surrounding biofluid are modeled using field 
theory. In the model, the biomatter is coupled to the gauge field representing the biofluid. 
It is shown that at non-relativistic limit various equation of motions, from the well-known 
Sine-Gordon equation to the simultaneous nonlinear equations, can be reproduced within 
a single framework. 
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1. Introduction 

The deoxyribonucleic acid (DNA) is getting the very important biomolecules. Es- 
pecially its helical structures undergoes a very complex dynamics and plays several 
important roles in various biological phenomena such as storage of information, in- 
heritance (replication, etc) and the usage of genetic information (transcription, etc) . 
Physically, transcription processes begin with the open state in the DNA double he- 
lix. The phenomena describe the internal DNA mobility such as rotational motion 
of nitrous bases, asymmetry of the helix, strength of the hydrogen bond and so on 
2} The rotational motion mode with such asymmetry has a solitary wave solution 
describing the open state of double helix ^ 

DNA is not motionless. It is in a constantly wriggling dynamics state in a 
medium of bio-organic fluid in the nucleus cell. One can think of the DNA molecule 
as behaving somewhat like a microscopic piece of worm-like elastic suspended in 
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solution. In this approximation, one ignores all details o f th e chemical structure. 
The motion of DNA surrounded by fluid is rarely studied Previous studies are 
usually done by solving the fluid equations and its wave equations simultaneously 
using appropriate boundary conditions H. For example, the thermal denaturation 
of double strand DNA depends on the solution surrounding the DNA molecules. It 
is necessary to take into consideration that the solving water does act as a viscous 
medium and may damp out bubble denaturation. The impact of viscosity has been 
investigated by Zdrakovic et al ^ which shows that the behavior of DNA dynamics 
in viscous solution is described by the damped nonlinear Schrodingcr (NLS) equa- 
tion. The solution is then obtained by expansion and performing order-by-order 
calculation. In these approaches, anyway the picture of interaction between DNA 
and its surrounding fluid is not clear. Also, in most models the over-damped DNA 
dynamics are treated by putting some additional terms by hand in the differential 
equation to obtain the non-homogeneous ones. 

There is a technology for DNA manipulation developed in the last decade using 
gel electrophoresis^. This technique is used to estimate the size of genetic sequences 
and understand the dynamics of individual strands of DNA as they move through 
a fluid solution. In theoretical polymer dynamics, the study is mostly focused on 
describing single-polymer dynamics to improve models of bulk properties. All of 
them attempt to understand the behavior of individual polymers under controlled 
fluid flows through the Navier-Stokes equation coupled with elastic strings. Theo- 
retical studies of this matter is still an open problem, especially related to reological 
problem. 

In this paper, we study the soliton excitation in DNA with surrounding fluid 
based on the gauge filed theory formulation El. The DNA dynamics is modeled as 
the result of interactions among matters in a fluid medium using the relativistic 
and gauge invariant fluid dynamics lagrangian. The theory is a relativistic one, 
but we can take its non-relativistic limit at the final stage to deal with problems 
in biomatter as done in some previous works, like in some models using the ideal 
gas approximation, Sine-Gordon models etM The la grangian is originally devoted 
to model the magnetofluid in plasma as a relativistic fluid system inspired by the 
similarity between the dynamical properties of fluid and electromagnetic field EI. 
Within the model, the interactions of DNA and fluid dynamics are described in a 
general way as the results of interactions among the fluid and matter fields. 

The paper is organized as follows. In the first chapter the model of interaction 
between biomolecule and fluid is given. It is then applied to describe the dynamics 
of worm-like DNA molecule with varying velocity. The paper is finally ended by a 
summary. 

2. The model 

Some mechanical models of DNA have been proposed over the years which are 
focused on different biological, physical and chemical processes in which DNA is 
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involved. One of the important model related to the identification of the unwinding 
of double helix in a " bubble" is called the torsional Y model The lagrangian 
describing the bubble of Y model is represent by a scalar (boson) field governed by 
the bosonic lagrangian ^J^l, 

Anattcr = (^$) f + V (<S>) , (1) 

where V(&) is the potential. For example in the typical $ 4 — theory, 

V{*) = ±m%&*-^\(&9>)* , (2) 

where m$ and A are the mass of matter and the nonlinear constant The her- 
mite conjugate is $t = (<&*) for a general complex field $ describing a rotational 
(torsional motion) of DNA. 

The interaction in the Y model with fluid medium is described in the gauge 
field theory framework. In this theory, the above bosonic lagrangian is imposed 
to be gauge invariant under local (in general non-Abelian) gauge transformation, 
U ee exp[-iT a 6 a (x)} »1- iT a 6 a {x) with 9 a < 1 T a 's are generators belong to 
a particular Lie group and satisfy certain commutation relation [T a ,T b ] — if abc T c 
with f abc is the anti-symmetric structure constant. The matter field is then trans- 
formed as $ $' ee exp[— iT a 9 a {x)] $, with T a are n x n matrices while $ 
is an n x 1 multiplet containing n elements. It is well-known that the symmetry 
in Eq. (fT]) is revealed by introducing gauge fields A a ^ which are transformed as 

U%' = U*- - g {d^9 a ) + f abc 8 b lJc, and replacing the derivative with the 
covariant one, 2? M = + igT a U^. The fluids is represent by the gauge fields that 
guarantees the invariant properties of the system. Then the total lagrangian with 
some additional terms to keep its gauge invariance is, 

£ ^matter ~t~ ^gaugc ~t~ ^int j (3) 

where, 

'-gauge — ^ J ^v° ; V*J 

C int = -gJ*U al1 + g 2 (&T a T b <I>) ^ 

The strength tensor is S° v = d^U^ - d v U£ + gf abc U b U^, while the 4-vector current 
is, 

J* = -i[(d^T a 4>~&T a (d^)] . (6) 

The coupling constant g then represents the interaction strength between gauge 
field and matter. We should note that, however the current conservation is realized 
by the covariant current d^J^ = with J* = -i [(<9 Al $)1"T a $ - &T a (d^)] . 

The gauge boson {/ M is interpreted as a "fluid field" with velocity u^, and takes 
the formic, 



^ = ([/ o a ,U a )EE^0, 



(7) 
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with, 

uJ = 7 a (l,-v B ) ) (8) 

where <f> is an auxiliary boson field, while 7 = (l — |v a | 2 ) . Here we adopt the 
natural unit, i.e. the light speed c = 1. We should remark that the auxiliary field <f) 
is introduced to keep correct dimension. 

For microfluid one usually works with a single fluid. From the gauge theory point 
of view, this is realized by the Abelian gauge lagrangian. The total lagrangian in 
this case becomes, 

C = (d„$*) + hm% - ^A ($*$) 2 + g 2 U^U^ $*$ 

~ {d^U v - d v U^) (cW - d v U») + igU» [(d^*) $ - $* (^$)] , (9) 

using Eqs. ©, ©~© and @. The strength tensor is given by S£ v = d^-dJI^. 

The non-relativistic limit can be obtained by performing a transformation t — > 
t = it and putting 7—^1 respectively. Imposing the variational principle of action, 
one reaches at the Euler-Lagrange equation Therefore, the equation of motion 
(EOM) for $ reads, 

(a 2 -m|-2 5 2 |v| 2 )$ + ^A$ 3 = 0. (10) 

for a real $ field. This result leads to the well-known nonlinear Klein- Gordon equa- 
tion. On the other hand, the EOM of is given by, 

^ + (v.V)v = - 5 F, (11) 

where v is fluid velocity in non-relativistic limit. F is total external force (see '-'-'for 
detail derivation). This is the model the underlying model throughout this paper. 



3. DNA dynamics in the surrounding fluid 

Behavior of the DNA interaction with surrounding fluid depends on the boundary 
condition of fluid phenomena. One should deal with the boundary value problem for 
the fluid equation and then solve the DNA dynamics through Eq. (fTUj) . For the sake 
of simplicity, however let us fix the fluid flow in order to avoid solving the boundary 
value problem. 



3.1. Constant Velocity 

First of all, consider v = Vq is a constant and write down Eq. (|10[) in the standard 



form, 

„4<I. + A<J> : ' = .). (121 



dt 2 dx 2 
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in term of coordinate (t,x), where fh% = m% + 2g 2 v$. It's well known that the 
solution of the equation can be represented by the Jacobi elliptic function 
For A > the solution is, 



&(x, t) = Am<s>sn [/3(x + xq — Ct),m] , 



(13) 



where /3 = a/Ato|/((1 — C 2 )(l + to)) , A — ^/2m/ (1 + to), C is the phase velocity 
and to is a modulus. In the limit to — > 1, one obtains the usual kink solitary waves 
as follow, 



&{x, t) = m% tanh 



Ato| 



2(1 -C 2 ) 



(x + Xq — Ct) 



(14) 



For A < the solution is, 

$(a;, t) = Affi^dn [/3(x + xq — Ct),m] , 



(15) 



where /3 = y/ — Ato|/((1 — C 2 ) (2 — to)) and A = y/2/(2 — to). In the limit of m 
1, one obtains the pulse solitary wave solution, 



$(x,£) = TO^v^sech 



— Am 



1-C 2 



%{x + xo-Ct) 



(16) 



Therefore, one has two kind of solutions, namely the kink solution and the 
pulse solitary wave solution. These solutions can be usually used to describe the 
propagation of a bubble in a torsional motion of the Y model. 



^-m%^-2g 2 v\x,t)^ + X^ A =0, (17) 



3.2. Velocity as a small perturbation 

In this subsection, let us assume that the DNA-fluid interactions lead to small per- 
turbation of the homogeneous solution. In this case a linear perturbation technique 
can be used. One can rewrite the EOM as follow, 

<9 2 $ 5 2 $ 
dt 2 dx 

According to the linear perturbation theory, the solution of Eq. p"?]) can be treated 
perturbatively, 

$(x,t) = $ (x,t) +e® 1 (x,t) , (18) 

where e is a small parameter, i.e. v 2 (x,t) ~ e. Substituting Eq. (IT8l) into Eq. ([T7|) 
leads to the following equation, 

??-^?-"^ + A<i> ^ ' (19) 
for the zeroth order and, 



dt 2 dx 2 



y 2 (x,t)$i=0, (20) 
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for the first order. Here V 2 (x,t) 



4 - 2g 2 v 2 (x,t) +3A<f>^). Eq. JT^ is an 



ordinary nonlinear Klein-Gordon equation, while Eq. (|20[) is the linear Klein-Gordon 
equation with coefficient variable. 

<f>o is a solution of Eq. plfl) which is the term of Jacobi elliptic functions (Eqs. 
(JT3J) ~ ([TB]) ) . The solution is the function of variable x' = x — Ct. To solve Eq. 
(|20p and by assuming that v{x,t) = v(x'), one needs to introduce new variables, 
x' — x — Ct and t' = t— Cx In the new coordinate, Eq. (|20|) reads, 



+ V 2 {x')$ 1 = , 



(21) 



dt 12 dx' 2 

This result yields that the solution has the form of, 

$(x',t') = <t){x')e- lut ' . (22) 

Then substituting this result into Eq. pip , one immediately gets the Sturm-Liouville 
equation, 

d 2 ^ 



dx' 2 



+ V 2 (x') 







(23) 



where V 2 (x') = V 2 (x) — u 2 . In principles using the WKBJ expansion technique, 
the solution of this equation can be obtained. One should notice that the solution 
is exist only when the 'potential' slowly varies with x 1 . 

First, write V(ex') — V(X) with X is the scale and e<l. One has, 

''' 2 " V 2 {X)cp. 



and the solution is 
Co 



<f>(X,e) 



V(X) 



=exp 



dX 2 



V(Od( 







(24) 



Ci 



V(X) 



=exp 



(25) 

One can also study another case with fluid velocity having the form of U(x) = 
sin(x). 



3.3. Stratified velocity 

One of the important experiments is the DNA (biopolymer) dynamics in shear flow 
5. The experiment shows that in a steady shear flow, the molecules continually fluc- 
tuate between nearly fully coiled and fully stretched conformation while undergoing 
end-over-end tumbling. The fluctuations between conformations are aperiodic and 
exhibited a power law roll-off at a high frequency, which is attributed to a tumbling 
instability driven by Brownian motion. 

In the model, the calculation of shear effect on the dynamics within the Y model 
is easy to describe. The description of such DNA bubble dynamics in shear flow is 
geometrically depicted in Fig. [TJ 
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*> X 



Fig. 1. Geometry of torsional motion of DNA in the stratified fluids. 



In the microfluid experiments, the fluid is usually treated one dimensionally and 
assumed to be incompressible, non-inertia, fully developed, unidirectional, Newto- 
nian and steady flow. Under these assumptions, the EOM for fluid flows is^, 







dP 

dx 



1 d 2 v 

v dz 2 



(26) 



where v and P are the fluid viscosity and pressure. If the pressure gradient occurs 
only at one edge of channel with its length is namely L x , one can solve the equation 
using boundary condition v(z = —L/2) = v(z = L/2) = 0. Then fluid velocity in 
the parallel plate is given by, 



v(z) 



1 AP 

2v L T 



L2 

4 



(27) 



Note that the velocity distribution is parabolic in z variable and called a Poiseuille 
flow. 

The EOM for the DNA dynamics is then given by, 



(28) 
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Fig. 2. The rotational motion of the DNA in the term of a pulse solitary wave as a function of x' 
with the coupling constants g = (solid line) , g = 0.1 (dashed line) and rotational motion with 
small oscillation fluid velocity (dot line) for a fixed parameter set (m$, v, C, A) = (1, 5, 0.5, 0.2). 

where U 2 (z) = m% + 2g 2 v 2 (z). By solving the boundary value problem of fluid 
velocity, the solution is given by Eqs. (fT4")l and (IT51) for A > and A < respectively. 

4. Summary 

A new model for DNA dynamics immersed in fluid system has been introduced. 
The model is constructed within the gauge field theory framework to describe the 
interactions between DNA matter and surrounding fluid from the first principle. 
The dynamics is then studied for various types of fluid velocity. 

The basic behavior of the bubble dynamics is described by a pulse solitary 
waves in Eq. ((TH|) . The open state dynamics, that is usually described by the soliton 
solution, in the present model is realized in the case of uniform fluid velocity as 
depicted in Fig. [2] The open bubble should correspond to the nonlinear excitation, 
and thus it is induced by the nonlinear dynamics of DNA double helix itself H| 
From the figure, one can conclude that the bubble amplitude is getting greater as 
the fluid flows at a constant velocity. On the other hand, the dotted line represents 
the effects of torsional motion on the DNA immersed in the oscillating fluid, i.e. 
v = sin(u;x). However, the oscillation in the bubble dynamics is relatively small 
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Fig. 3. The rotational motion of DNA in term of a pulse solitary wave in the Poiseuille flow with 
a fixed parameter set (m$, v, C, A) = (1, 5, 0.5, 0.2) 

since by assumption the perturbation term (fluid velocity) is kept small. 

Another effects due to stratified flow in the bubble dynamics is shown in Fig. [3] 
As mentioned above, according to the experiments the average extension of DNA 
under shear flow was found to increase gradually with increasing velocity gradient, 
reaching a plateau below half-full extension. The hydrodynamics friction is more 
important factor to aligned of the DNA molecule in a direction parallel to the flow. 
The present result is given in Fig. [3] which shows that the fluid flow increases the 
soliton amplitude with the same profile, that is it reaches the peak at the center 
and almost unaffected near the plateau. 

Further, the effects of external forces like electromagnetic force might be inter- 
esting to b e investigated. This might be important as one considers the nanofluidic 
system H^. This will be worked out and discussed in the next-coming works. 

Acknowledgments 

AS thanks the DAP Consultant for partial funding and the Group for Theoretical 
and Computational Physics LIPI for warm hospitality during the work. This work 
is funded by the Indonesia Ministry of Research and Technology and the Riset 
Kompetitif LIPI in fiscal year 2011 under Contract no. 11.04/SK/KPPI/H/2011. 

References 

1. L. V. Yakushevic, Nonlinear physics of DNA 2nd edn (Wiley and Sons, 2002). 

2. L. V. Yakushevich, A. V. Savin and L. I. Mawvich, Phys. Rev. £66, p. 016614 (2002). 

3. D. Voet, J. G. Voet and C. W. Pratt, Fundamental of Biochemistry (John Wiley and 
Son, 2008). 

4. S. Zdrakovici, J. A. Tuszynski and M. V. Sataric, J. Comput. Theor. Nanosci, 21, 1 
(2005). 



4, 2012 2:43 WSPC/INSTRUCTION FILE sulaiman 



10 A. Sulaiman and L.T Handoko 

5. J. P. Rickgauer and D. E. Smith, Single-Molecule Studies of DNA (Kluwer Academic 
Pub., 2006), ch. Soft Matter: Scattering, Imaging and Manipulation. 

6. A. Sulaiman and L. T. Handoko, J. Comp. Theor. Nanosci. 8, 124 (2011). 

7. A. Sulaiman, A. Fajaruddin, T. P. Djun and L. T. Handoko, Int. J. Mod. Phys. A, 
24, 3630 (2009). 

8. L. H. Ryder, Quantum Field Theory 2nd Ed. (Cambridge University Press, 1996). 

9. F. Cooper, A. Khare, B. Mihaila and A. Saxena, Phys. Rev. E 72, p. 036605 (2005). 

10. M. W. Dingemans, Water Waves Propagation over Uneven Bottom (World Scientific, 
1997). 

11. M. Cadoni, R. D. Leo and S. Demelio, Int. J. Non-linear Mechanics 43, p. 1094 (2009). 

12. J. Viovy, Rev. Mod. Phys. 72, 813 (2000). 



